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Abstract 

A fit to proton, neutron and deuteron spin asymmetries is presented 
and polarized parton distributions in nucleon are given. These densi- 
ties have their roots in the MRS fit for unpolarized case. The integrals 
of polarized distributions are compared with the experimental figures. 
The role of polarized gluons is also discussed. 
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There are a lot of new experiments which yield the data on spin structure 
of a nucleon. We have the newest experimental points from CERN [|l], |^ and 
SLAC ^, and also the older ones from these laboratories [Q, [0. The 
spin asymmetries are measured on proton (SLAC- Yale 0, EMC 0, SMC 
[|] and E143 @]), neutron (E142 g made on ^He) and deuterium (SMC 
[H] and E143 0) targets. Using the experimental figures we can study phe- 
nomenologically nucleon spin structure and in particular we can determine 
the polarized parton distributions. 

For unpolarized case several fits were performed [8-11], there were also 
attempts to get the spin distributions [12-18]. Some time ago Martin, Roberts 
and Stirling (MRS) [g presented a complete fit with determination of parton 
(i.e. quark and gluon) distributions. Quite recently after the measurements 
from Hera [|l^ for the small x region the new fits with more reliable gluon 
contribution and modified sea were performed [Q. Our discussion how to 
determine the polarized structure functions using MRS fits |0] was given in 
Ref.fli and in Ref.^. 



In this paper we would like to get polarized quark parton distributions 
starting from the unpolarized ones and using existing data on spin asymme- 
tries. The calculated values of octet axial-vector couplings 03 and as obtained 
from the fit are compared with the experimental values gotten from nucleon 
and hyperon /3-decays and modified for QCD corrections. We use these quan- 
tities to find the best fit (of course together with the values). In order to 
stabihze the fits we put the restriction on as (such procedure takes place also 
e.g. in Ref . fT^ ) . We have also tried to include polarized gluons contribut- 



ing in the way proposed in Ref.[T^. It comes out that the polarized gluonic 
degrees of freedom do not lead to any substantial improvement in a fit. 

Let us start with the formulas for unpolarized quark parton distributions 
given (at = 4 GeV^) in the newest fit performed by Martin, Roberts and 
Stirling ||^. We have for the valence quarks: 

u^{x) = 2.704a;-°-^°^(l-a;)^-^'^(l-0.76v^ + 4.20x), 

d^{x) = 0.251x-°-^''^(l-x)^-^^(l + 8.63v^+0.32x), (1) 

and for the antiquarks from the sea: 

2u{x) = 0.392M(x) -(5(x), 
2d{x) = 0.392M{x) + 5{x), 
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2s(x) = 0.196M(x), 
2c{x) = 0.020M(x). 



(2) 



In eq.(2) the singlet contribution is: 

M{x) = 1.74x-i-°^^(l - xy°-\l - 3A5Vx + 10.3a;), (3) 
whereas the iso vector part: 

S{x) = 0.043x~°-^(l - xY^-\l + 64.9x). (4) 
For the unpolarized gluon distribution we get: 

G{x) = 1.51x-^-3°^(l - xf-^%1 - A.U^ + 10. Ix). (5) 

We assume, in an analogy to the unpolarized case, that the polarized 
quark distributions are of the form: x"(l — x)^P2{y/x), where P2{y/x) is 
a second order polynomial in ^/x and the asymptotic behaviour for x^O 
and X— s>l (i.e. the values of a and (3) are the same (except for AM, see a 
discussion below) as in unpolarized case. Our idea is to split the numerical 
constants (coefficients of P2 polynomial) in eqs.(l, 3 and 4) in two parts in 
such a manner that the distributions are positive defined. Our expressions 
for Aq{x) = q~^{x) — q~{x) (g(x) = q'^{x) + q~{x)) are: 

Au^{x) = x-^-^°'^{l-xY-^%ai + a2Vx + a3x), 

Ad^{x) = x~^-^^\l-xf-^\bi + b2V^ + b3x), 

AM{x) = x-°-^^'(l-x)i°-i(ci + C2v^), (6) 

A5{x) = x-°-^(l-x)^°-^C3(l + 64.9a;). 

For a moment we will not take into account polarized gluons i.e. we put 
AG = 0. For total sea polarization i.e. AM, we assume that there is no 
term behaving like x"^"^^ at small x (we assume that AM and hence all 
sea distributions are integrable), which means that coefficient in this case 
have to be splitted into equal parts in M"*" and M~ . The next term (x~°'^^^) 
is relatively more singular then the sea contribution in the preferred fit in 
Ref. ||17||. On the other hand when we put the coefficient in front of x""'^^^ 
equal to zero (i.e. Ci = 0) M{x) will behave for x — like x"*^ *^^^. In this 
case we get the fit with higher x^/^df- That means that in spite of the fact 
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that in the present fit the unpolarized sea behaviour is less singular for small 
X values, contrary to the case in Ref . [ ITj] , the model for AM{x) with more 
singular sea contribution is phenomenologicaly chosen. As we will see that 
influences the behaviour of gi{x) for small x values. 

In order to get the unknown parameters in the expressions for polarized 
quark distributions at = 4GeV^ (see eq.(6)) we make a fit to the experi- 
mental data on spin asymmetries for proton, neutron and deuteron targets. 
The theoretical expressions for these asymmetries are given by: 

_ 4Au^{x) + Ad^{x) + 2.236AM(x) - 3A6{x) 
'^^^ ~ 4u,{x) + d,{x) + 2.236M(x) - 36{x) ^ ^ 
,n(. ^ An,(x)+4A4(x)+2.236AM(x)+3A^(x) 

u,{x)+Ad,{x)+2.23QM{x)+36{x) ^ ^' ^ 

^ 5Au,{x) + 5Ad, + 4.472AM(a:) 3 ^ p. 

= 5u.{x) + 5d. + 4A72M{x) ' + R). 

Where the ratio R = a^/ ax (which vanishes in the Bjorken limit) is taken 
from and pn'm a. probability of D-state in deuteron wave function (equal 
to (5±1)% §§. 

We assume that the spin asymmetries do not depend on Q"^ (it is only our 
first order approximation) what is suggested by the experimental data [|T], Q 
and phenomenological analysis . We hope that numerically our results at 



= 4 GeV^ will not change much if the evolution of Fi and gi functions 
will be taken into account. 

Spin structure function e.g. is given by: 

gl{x) = (4Am^(x) + Ad^{x) + 2.236AM(x) - 3A5(x))/18. (8) 

The obtained polarized quark distributions Au{x), Ad{x), AM{x) and 
A6{x) can be used to calculate first moments. For a given we can write 
the relations: 

4 , 1 , , 1 , 4 , 

r? = —Am H Ad H As H Ac, 

^ 18 18 18 18 

14 14 
r" = —An + —Ad + — As + —Ac, (9) 
^ 18 18 18 18 ' ^ ^ 

where Aq = Aq{x) dx and Fi = gi{x) dx. 
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We define other combinations of integrated quark polarizations: 



03 = Au — Ad, 

as = Au + Ad-2As, (10) 
AS = Au + Ad + As, 

Sucli results for the integrated quantities (calculated at 4 GeV^) after 
taking into account known QCD corrections (see e.g. Ref.|Q) could be com- 
pared with axial-vector coupling constants Qa and Qs known from neutron 
/3-decay and hyperon /3-decays (in the last case one needs SU{3) symmetry). 
The difference of rf((5^) and T^IQ"^) can be expressed by: 

r?(g')-r^(Q') = {Aum-Ad{Q'))/6 = a^m/G = cns{Q^)9a/Q, (H) 



where cjqsiQ"^) describes QCD corrections for non-singlet quantities and 



qa = 1.2573 ±0.0028 (see e.g. Ref . P3| ) is obtained from the neutron /3-decay. 
In our paper is constant and takes the value 4GeV^. 

We get from the experimental figure 03(4 GeV^) = cns{.'^ GeY'^)gA = l-H 
and with this value we shall compare calculated from our fits. Another 
useful combination of V^i^Q"^) and V'^{Q'^) is equal to: 

r?(Q2) + r^(Q2) ^ 5a8(Q2)/i8 + 2As(Q2)/3 (12) 

with og = ciys{Q'^)g8, where gs = 0.58±0.03 p3| is obtained from the hyperon 
/3-decays. Knowing c^siQ"^) we can calculate 05(4 GeV^) = 0.51 ± 0.03 and 
with this number we shall compare the results obtained from our fit. If we 
have had very precise experimental data and in the whole x range there 
would be no problems with determination of polarized quark distributions. 
Unfortunately that is not the case yet. Actually, from the experiment we have 
information on Fj* and F". The combination Fj'-F" is directly connected to 
qa experimental quantity modified by QCD corrections. On the other hand 
Fj'+F" is the combination of ag and As and it came out in Ref. |r?[] that the 
fits are not sensitive enough to determine and As separately in a stable 
way. The values of ag and As were different for our models and for different 
subsets of data. To stabilize the determination of parameters we also here 
assume in addition that ag = 0.51 with 0.1 as artificial theoretical error. 
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We get the following values of our parameters (describing the polarized 
quark distributions in eq.(6)) from the fit to all existing data for spin asym- 
metries: 

ai = 1.71, 02 = —6.68, = 14.2, 

= -0.005, 62 = 0.835, h = -3.33, (13) 
ci = -1.10, C2 = 1.38, C3 = -0.03. 

In the second row of the Table 1 the integrated quantities: Ff, F", 03, ag, 
AS and AM together with the coresponding x^/^df that follow from our 
fit are presented. They can be compared (first row) with the values obtained 
from our previous fit presented in Ref. ||171|. We will not show the comparison 
of the results of the fit to CERN data only (proton + deuteron) and SLAG 
data (proton + neutron + deuteron). Two fits like before are in agreement 
with each other and with the overall fit. It is specially interesting because 
the CERN data are taken at much smaller x values then the SLAG data. 

In Figs.(l, 2 and 3) we present the comparison of our fit with the exper- 
imental asymmetries for proton (1), neutron (2) and deuterium (3) target. 
We see that in the case of deuterium many experimental points lie below our 
fitted curve especially in the large x region where the errors are relatively big. 
In the Fig. (2) we see that A" approaches very slowly zero for x ^ what 
produces relatively high (negative) value of F". Because we have 62 points 
for proton spin asymmetries in comparison with 33 for deuteron and 8 for 
neutron it seems that our curves are dominated by the proton data. Small 
discrepancies between different experiments are also not excluded. 

The obtained quark distributions lead to the following integrated quanti- 
ties: Am = 0.69 {Au^ = 0.83), Ad = -0.45 (A^ = -0.08) and As = -0.13, 
hence the amount of sea polarization is AM = —0.65. This last number is 
not small and hence also the strange sea polarization in our model is rather 
big. As one can see from Table 1 the 03 value seems to be in our fit very 
close to the predicted value 1.11. AS is rather small and the tendency that 
the model with more singular sea behaviour (ci 7^ 0) produces very small 
AS values is confirmed (this was already observed in Ref.|T^). Our inte- 
grated quantities (F^ and F") differ slightly from the values quoted by the 
experimental groups, whose figures are calculated directly from the experi- 
mental points with the assumption of Regge type behaviour at small x. But 
on the other hand our polarized quark distributions satisfy all the constraints 
taken implicitly into account in fits to the unpolarized data. The difference in 



Table 1 
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integrated quantities comes mainly due to our assumptions about small x be- 
haviour for spin-dependent distributions. If we compare our results integrated 
in the interval 0.029 < x < 0.8 (region covered in the SLAG experiments) we 
get rj'=0.11 (SLAG result without the extrapolation to unmeasured region 
is 0.12±0.01) and Tf = 0.04 (to be compared with 0.04±0.005). Also we use 



R pO| which is not very realiable for small x and this also gives an additional 
error, which is rather difficult to estimate. 

It is interesting to see what will happen when we relax our condition 
as = 0.51 ± 0.1. We get the fit with comparable P^r degree of freedom 
(x^/Ndf = 1.02) and = 0.12, = -0.08, ag = 1.23, as = 0.23, AS = 
0.13 and AM = —0.18. We see that ag value ( Og ~ AE ) is rather small in 
this case and also we have much smaller sea polarization. 



We have also tried to include polarized gluons along the line of ref. JTS 
assuming for the gluon distribution: 

AG{x) = x-°-^°i(l - xf^^di + rf2v^), (14) 

with a new di and d2 constants which have to be fitted. The appearance 
of non-zero gluonic distribution affects our formulas only through the sub- 
stitution: Ag ^ Ag - f^AG. In such fit {x^/Ndf = 1-01) we got (after 
integration) the negative sign of the gluonic contribution (if we take gluonic 
distribution with the full strength at x = 0). In this case Ff is equal to 0.19 
because the gluonic contribution is added to the quark contribution instead 
of being substracted. 

In the row three of Table 1 the results of the fit with gluons are shown 
in which gluon contribution is less singular AG{x) ~ for small x i.e. 

the coefficient in front of is equaly devided between G^{x) and G^{x) 

{di = 0). In this case I^df = 1-02 and is comparable to the fit without 
gluons. In this case sea polarization is very small and also AE, as expected, 
is rather small (0.06). Also we get in this case for the value higher than the 
experimental figure. We conclude that the inclusion of gluonic contribution 
does not lead to the substantial improvement of the fit. 

The Fig. (4) shows the comparison of gi{x) calculated from our fit with 
the experimental points (evolved to common value = 4GeV^). We do 
not observe the growth of Qi for small x values in our model contrary to the 



previous fit in Ref. [17|. It is caused by the relatively singular behaviour of 



the sea contribution for small x values. 
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Starting from the new, improved version of the MRS fit ||^ to the un- 
polarized deep inelastic data we have made a fit to proton, neutron i^He) 
and deuteron spin asymmetries in order to obtain polarized quark parton 
distributions. To stabilize the fits we added the experimental information on 
octet quantity Og. We have calculated the parameters of the polarized quark 
distributions using the combination of all existing proton, neutron and deu- 
tron spin asymmetries measurements (including the new results on proton 
and deutron E143 experiments and very recent SMC deuteron data from 
CERN). We do not need gluonic contributions to be taken into account, i.e. 
the fit with gluons is not better. The next step in front of us is to include Q"^ 
dependence of spin asymmetries in comparison with experimental data. 
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Figure captions 



Figure 1 The comparison of spin asymmetry on protons with the curve gotten 
from our fit to all existing data. Points are taken from SLAC (E80, 
E130, E143) and CERN (EMC, SMC) experiments. 

Figure 2 The comparison of spin asymmetry on neutrons (SLAC E142 data) 
with the curve gotten from our fit. 

Figure 3 Our prediction for deuteron asymmetry compared with the SMC and 
SLAC data. 

Figure 4 The data for gi{x) structure function with the curve gotten using the 
parameters of our fit. The data points are taken from SLAC and CERN 
and evolved to common Q'^ — 4GeV^. 
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Table 1 



The first moments of polarized distributions (see eqs.(9) and (10)). The 
strange sea polarization As (not presented in the Table) is connected to 
the total sea polarization by the relation: As = 0.196AM. We have made 
our fits taking all existing experimental data on spin asymmetries. 



Fit 


r? 


1 


0.3 


as 


AS 


AM 


xVNdf 


Old fit [17] 


0.139 


-0.072 


1.27 


0.47 


0.20 


-0.45 


1.05 


New fit 


0.119 


-0.072 


1.14 


0.50 


0.12 


-0.65 


1.01 


New fit with gluons 


0.123 


-0.082 


1.23 


0.51 


0.06 


-0.09 


1.02 
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